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Abstract. Spherical Hartree-Fock calculations with projection onto zero total linear momentum before 
the variation are performed for the nuclei "‘He, ‘^C, ‘®0, ^®Si, and "‘°Ca using a density-independent 
effective nucleon-nucleon interaction. The results are compared to those of usual spherical Hartree-Fock 
calculations subtracting the kinetic energy of the center of mass motion either before or after the variation 
and to the results obtained analytically with oscillator occupations. Total energies, hole-energies, elastic 
charge form factors and charge densities and the mathematical Coulomb sum rules are discussed. 

PACS. 2 1.60.-n Nuclear-structure models and methods 


1 Introduction 

We consider the nucleus as a closed system of interact¬ 
ing, non-relativistic nucleons. The homogenity of space 
requires that the total linear momentum of this system is 
conserved. Consequently the hamiltonian describing any 
particular nucleus cannot depend on the center of mass 
(COM) coordinate of its constituents, but (besides on spin- 
and isospin-quantum numbers) only on relative coordi¬ 
nates and momenta. The dependence on the total momen¬ 
tum is trivial: it describes the free motion of the total sys¬ 
tem and can always be transformed away by considering 
the system in its COM rest frame. We have then to solve 
the corresponding Schrodinger equation for the remaining 
“internal” hamiltonian. In principle this can be achieved 
by writing this hamiltonian in .Jacobi coordinates. How¬ 
ever, nucleons are fermions and thus do obey the Pauli 
principle. Since the Jacobi coordinates depend on all the 
nucleon coordinates, thus an explicit antisymmetrization 
of the wave functions is required as it is performed, e.g., in 
few-body physics. Being already there sometimes rather 
involved though still feasible, such an explicit antisym¬ 
metrization becomes impossible in the many-body system 
(e.g., the antisymmetrization of 20 like nucleons would 
require 20 factorial different terms). Thus in the many- 
body system the antisymmetrization usually is performed 
implicitely by expanding the wave functions in terms of 
Slater (or generalized Slater) determinants. In this way 
the Pauli principle is automatically fulfilled. The Slater 
determinants, however, depend on 3A instead of the al¬ 
lowed 3A-3 coordinates and thus contain contaminations 
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due to the motion of the system as a whole, so called “spu¬ 
rious” admixtures. Galilei-invariance is broken. 

This defect of almost all microscopic nuclear struc¬ 
ture models has been recognised almost immediately 
after the development of the shell-model. It was shown 
later on [2 that in case of pure harmonic oscillator con¬ 
figurations one can get rid of this problem by diagonaliz¬ 
ing the (oscillator) COM hamiltonian and projecting all 
states not corresponding to the ground state of this op¬ 
erator out of the spectrum of the many nucleon hamilto¬ 
nian. This procedure, however, requires the use of so called 
complete nfiw-spaces (since only then COM and internal 
excitations decouple exactly) and thus is of little help in 
most of the usual approaches to the nuclear many-body 
problem. A more general solution is the projection of the 
wave functions into the COM rest frame jH], which ensures 
translational- and, if performed before solving the cor¬ 
responding Schrodinger equation (usually by variational 
methods) even full Galilei-invariance 0. The key idea of 
this projection is to superpose the wave function shifted 
all over normal space with identical weights and thus to 
achieve vanishing total linear momentum. Since the bound 
states of a nucleus are localized, this procedure always 
does converge (for scattering states a slightly different pro¬ 
cedure has to be used USD- The projection method has the 
advantage that it works in general model spaces as well as 
for general (non-oscillator) wave functions. 

Though in principle known since almost half a century, 
only few practical calculations have been performed using 
this method. The reason for this is quite simple : the pro¬ 
jection operator is an A-body integral operator with the 
rather nasty property to link the usual model space states 
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Fig. 1. In the left part the total binding energy of ^He as obtained with the Brink-Boeker force B1 [14] using an oscillator length 
b=1.41 Fm is plotted against the size of the basis. Three curves are shown : the open circles correspond to a normal spherical 
Hartree-Fock calculation with the expectation value of /2MA subtracted after convergence, the full triangles display the 
results if this center of mass correction is included during the iteration, and, finally, the full circles give the result of a spherical 
Hartree-Fock calculation with projection into the momentum rest frame before the variation (CMPSHF). Furthermore, the figure 
displays the results of constrained calculations adding the oscillator center of mass Hamitonian with a Lagrangian multiplier A 
to the internal hamiltonian which “penalizes” COM excitations. The right part of the figure displays the hole energies obtained 
in the oscillator limit HO (Nmax=0) and in the Hartree-Fock approximation HF (with Ninax=18) using the “normal” (no) as 
well as the “projected” (pr) approach in both cases. Note that the normal approach includes the usual COM correction (see 
text). 
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to rather highly excited (and thus usually unoccupied) 
ones as well as to the fully occupied ones, which are often 
treated as an inert core. This is easy to understand : any 
change of the linear momentum of the valence nucleons 
requires a corresponding change of the linear momentum 
of the core in order to ensure vanishing total linear mo¬ 
mentum for the system. Unlike angular momentum, linear 
momentum is thus a true A-body correlation and hence 
much more complicated to treat than the latter. 

Because of these difficulties, instead of treating Galilei- 
invariance correctly, its breaking is usually neglected adopt¬ 
ing the well known text book argument that it induces 
only 1/A effects and thus “can savely be neglected for nu¬ 
clei beyond oxygen” 0, provided the usual approximate 
corrections like subtracting the kinetic energy of the COM 
motion from the original hamiltonian or the use of the so 
called Tassie-Barker factor in the analysis of form fac¬ 
tors are done. 

That this, however, is not true has been shown by sev¬ 
eral studies within the last decade. Hartree-Fock calcula¬ 
tions with projection into the COM rest frame for "^He [6] 
as well as the analysis of form factors and charge densities 
of several spherical nuclei [711^ have demonstrated that 
the correct treatment of Galilei-invariance yields consider¬ 


able effects far beyond the usually assumed 1/A level. The 
same holds for scattering states as demonstrated in ref. 
^ for the inclusive quasi-elastic electron scattering again 
from ^He. Recently now, a whole series of model investiga¬ 
tions |1 Din 11112111 .^l has been published, in which the COM 
effects have been studied in a more systematic way. Con¬ 
siderable effects have been seen for spectral functions and 
spectroscopic factors, transition form factors and densi¬ 
ties, energies of hole-states. Coulomb sum rules, response 
functions and many more. These investigations, however, 
have been undertaken with rather simple wave functions 
: the ground states of the doubly even A-nucleon systems 
^He, ^®0 and “^^Ca have been described in the simple os¬ 
cillator limit and for the ground and excited states of the 
corresponding odd (A-l)-nucleon systems simple one-hole 
states have been used. This has the advantage that all cal¬ 
culations can be performed analytically but is definitely 
not very realistic. So, e.g., the above mentioned pure oscil¬ 
lator A-nucleon configurations are non-spurious and thus 
the projection yields here no additional effect with respect 
to the usual approach to subtract the kinetic energy of the 
COM motion. It is hence desirable to study, e.g., these 
ground states in more realistic approaches. This will be 
done in this and a forthcoming paper. 
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2. Same as in Fig. 1, but for the nucleus with oscillator length b=1.72 
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Fig. 3. Same as in Fig. 1, but for the nucleus with oscillator length b=1.79 Fm. 
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Fig. 4. Same as in Fig. 1, but for the nucleus ^®Si with oscillator length b=1.85 Fm. 


For this purpose we have performed spherical Hartree- 
Fock calculations with projection into the COM rest frame 
before the variation for the six nuclei "^He, ^®Si, 

and ^°Ca. The results have been compared with those 
of normal spherical Hartree-Fock calculations subtract¬ 
ing the kinetic energy of the COM motion either before 
or after the variation and with the analytically obtained 
oscillator results out of ref. m For each of the consid¬ 
ered nuclei up to 19 major oscillator shells have been used 
as single particle basis. As effective interaction the sim¬ 
ple Brink-Boeker force B1 m has been taken. We are 
aware of the fact that this interaction is not very realis¬ 
tic. However, the aim of the present investigation is not 
a comparison with experiment but the study of the ef¬ 
fects of a correct treatment of Galilei-invariance. For this 
purpose the B 1-interaction is as good as any other. Fur¬ 
thermore, consisting out of Gaussians, it can be treated 
in the oscillator limit analytically and thus allows for a 
direct comparison with the results reported in ref. m 

Section 2 of the present paper gives a short summary of 
the spherical Hartree-Fock approach with projection into 
the COM rest frame before the variation. Section 3 will 
then describe some details of the calculations and present 
the results for the total energies, the hole-energies, the 
elastic charge form factors and corresponding charge den¬ 
sities and the Coulomb sum rules. Conclusions, three ap¬ 
pendices with some detailed formulas and references con¬ 
clude the present paper. 

In the second of the present series of two papers we 
shall then discuss the effects of the correct treatment of 
Galilei-invariance on the spectral functions and spectro¬ 
scopic factors obtained with the wave functions out of the 
present paper. 


2 COM-projected Hartree-Fock. 

The essential mathematics for Hartree-Fock calculations 
with projection into the COM rest frame before the vari¬ 
ation has been presented in detail already in ref. [6] and 
hence will be summerized only briefly in the following. We 
start by defining our model space by Mf, oscillator single 
particle states, the creators of which will be denoted by 
{cl; i = 1, ..., Mb}. We shall furthermore assume that the 
effective hamiltonian appropriate for this model space is 
known and can be written in the chosen representation as 
a sum of only one- and two-body parts 

H = '^t{ir)clcr -t- '^v{ikrs)clclcsCr (1) 

ir ikrs 

where t{ir) are the single particle matrix elements of the 
kinetic energy operator and v{ikrs) the antisymmetrized 
two-body matrix elements of the considered interaction. 
We shall assume that this interaction is translational in¬ 
variant, i.e., it does not depend on the center of mass 
coordinate of the two nucleons. Density dependent inter¬ 
actions (in their usual form) do not fulfill this requirement. 
Their treatment is much more complicated as has been de¬ 
scribed in detail in ref. m Such interactions will not be 
considered in the present paper. 

In the Hartree-Fock approach one searches for the op¬ 
timal one-determinant representation of the A-nucleon 
ground state having the form 

i^) = |n^i}io) (2) 
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5. Same as in Fig. 1, but for the nucleus ®^S with oscillator length b=1.91 Fm. 


where 


and 


0/3 


E 


i=l 


^i/3‘ 


t 

i 


of the total momentum of the considered system) with 
identical weights. 

(3) The energy functional Epr to be used in the Hartree- 
Fock approach with projection into the COM rest frame 
before the variation can then be written as 


Mb 

— 'y ^ DipCi (4) 

i=l 

respectively, with D being a unitary (M/, x M/,) transfor¬ 
mation. In eq. (2) we have assumed that the selfconsis- 
tent states created by the operators (3) are ordered ac¬ 
cording to their energy so that (3 = 1,..., A correspond to 
the occupied “hole” states h, h', .... The unoccupied or¬ 
bits f3 = A+ 1,..., Mb will be denoted as “particle” states 
p, p', ... in the following. 

Now, obviously, the determinant 0 is not translation- 
ally invariant. In order to obtain a Galilei-invariant wave 
function we have to use instead of m the expression 

IZ,; 0, = (5, 

^/{D\C(0)\D} 

as test wave function in the variation. Here 

(7(0) = jSaS{a) (6) 

with 

S{a) = exp{/a • P} (7) 


_ /da a^fdng h(a) 
f da f dfia n{a ) 

where we have introduced the shifted overlap function 

n{a) = {D\S{a)\D) = det7f(a) (9) 

which can be represented as the determinant of an (Ax A)- 
matrix 

Xhh'ia) = {h\S{a)\h') = ^ DihSik{a)Dlb, (10) 

ik 

with Sik {a ) being the matrix representation of the shift 
operator within the chosen harmonic oscillator single par¬ 
ticle basis. These matrix elements are given in the ap¬ 
pendix A. Furthermore we use in 0 the shifted energy 
function 

h{a) = {D\HS{a)\D) = t{a) + V 2 {a) (11) 

Here the one-body term is given by 

t{a) = n{a) ^ t{ir) pri{a) (12) 

ir 

with the shifted density matrix being defined as 


projects into the COM rest frame by superposing all states 
created by the shift operator 0 (here P is the operator 


^ Srkia) ^ Dlb [X-\a)U' A/z' (13) 

k hh' 
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Fig. 6. Same as in Fig. 1, but for the nucleus ^®Ca with oscillator length b=1.97 Fm. 
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Fig. 7. In the left part of the figure the square of the charge form factor for ^He is displayed as function of the 3-momentum 
transfer. For the nucleon form factors the usual dipole parametrisation [21] has been used. Open inverted triangles correspond to 
an oscillator occupation with no COM correction included (normal oscillator), full triangles give the oscillator result including 
the Tassie-Barker factor [5] (projected oscillator), open circles display the form factor obtained with normal Hartree-Fock 
(including the COM correction in the hamiltonian during the variation) taking into account the dynamic correction (essentially 
again the Tassie-Barker factor) out of the text. Finally, the full circles display the result of the projected calculation. The right 
part of the figure gives the corresponding charge densities (obviously calculated for point nucleons). 
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Fig. 8. Same as in Fig. 7, but for the nucleus 


and for the two-body part of m one obtains 


for the occupied and 


V 2 {a) = n{a) v{ikrs) psk{a) Priia) (14) 


ikrs 


The energy functional m has to be minimized with 
respect to arbitrary variations of the underlying Hartree- 
Fock transformation D. This transformation, however, has 
to be unitary and thus not all of the Mb ■ Mb matrix el¬ 
ements of D are linear independent. Nevertheless an un¬ 
constrained minimization of the functional (8) can still be 
performed, if one parametrizes the underlying Hartree- 
Fock transformation D via Thouless’ theorem m, which 
states that any Hartree-Fock determinant \Dd) can be 
represented in terms of the creation and annihilation op¬ 
erators of some reference determinant \Do) via 


iDd) 


c{d) exp 



dphbp{Do)bhiDo) 


\Do) 


(15) 


provided that the two determinants are non-orthogonal, 
since 

c(d) = {Do\Dd) (16) 


bliDd) = - E dph'bliDo) 

p' V h' 

(18) 

for the unoccupied states, respectively. They are given in 
terms of the {Mb -A)-A linear independent variables dpt ■ 
The {A X A) matrix L in 111 711 is defined by the expression 

1a + d^d* = W (19) 

while the {{Mb — xl) x {Mb — 7l)) matrix M out of 111 8ll 
can be obtained by the solution of the equation 

1m,-a + d*d^ = MMl (20) 


The variational equations resulting from the minimization 
of the functional m thus get finally the form 


dEpr 

ddpb, 




= 0 


( 21 ) 


where the {{Mb — A) x A) matrix G is defined as 


Gph = 


J da J dGg gph{a) 
J da a'^ J dQa n{a) 


( 22 ) 


and the function gph{a) is given by 


The creation operators belonging to the Hartree-Fock de¬ 
terminant \Dd) are then related to those of the reference 
determinant |Zlo) via 


bliDd) = 

h' 



E 



gph{a) = {D\ [H - Epr] S{a)bl{D)bh{D) \D) = 

Mb A f 

[X~^{a)]hh'Dih'<[h{a) - Eprn{a)] Sir + 
i,r=i h' = l I 

Mb 'j Mb 

”(“)E Sr, {a) Dip 

k=l } s=l 


(17) 


(23) 
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Fig. 9. Same as in Fig. 7, but for the nucleus 
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where 

Psr (a) (24) 

rs 

The “local” gradient vector ll^ . obviously, has to vanish 
at the solution of CB, too. This solution can be obtained 
using standard methods as they have been described, e.g., 
in Ref. [T^ . 

Up to now no symmetry restrictions have been im¬ 
posed on the Hartree-Fock transformation Q, 0- Thus 
the Hartree-Fock vacuum 0 breaks in general besides the 
translational invariance also other symmetries like, e.g., 
the conservation of the total angular momentum and the 
parity. For this general case therefore besides the momen¬ 
tum projection also the projection on these other symme¬ 
tries would be required. The situation becomes, however, 
much simpler, if only spherically symmetric Hartree-Fock 
transformations are admitted. Then each of the selfcon- 
sistent states created by the operators m has the isospin 
3-projection, the orbital and total angular momentum and 
the 3-projection of the latter as “good” quantum numbers, 
and the sums in m and 0 run only over the node quan¬ 
tum number. For nuclei with closed angular momentum 
subshells the corresponding Hartree-Fock vacuum 0 has 
then total angular momentum = 0+ and conserves the 
proton as well as the neutron number. Consequently the 
projection on these symmetries becomes redundant and 
we are left with only the linear momentum projection as 
described above. 

Furthermore, for spherically symmetric systems obvi¬ 
ously neither the shifted overlap 0 nor the correspond¬ 
ing energy function 111 111 do depend on the direction of the 
shift vector a. Thus the angle integrations in © and CIJ 
induced by the operator 0 become trivial and only a sin¬ 
gle integral over the radial component of the shift vector 


remains to be done numerically. An explicit formulation 
of this special case will not be given in the present paper. 
However, it is obvious that the calculation of the expres¬ 
sions needed for the minimization of the energy functional 
® is then simplified considerably. 


3 Results and discussion. 

We have considered the six nuclei "^He, ^®Si, 

and ^°Ca. As hamiltonian, as in ref. ca, the Brink- 
Boeker interaction B1 m complemented with a short 
range (0.5 Fm) two-body spin-orbit term having the same 
volume-integral as the corresponding zero-range term of 
the Gogny-force DIS na, plus the Coulomb force and 
the kinetic energy has been used. First, the energy of 
the simple oscillator determinants for these nuclei (e.g., 
(Os)^(Op)^^ for ^®0) has been minimized with respect to 
the oscillator length-parameter b. For the intermediate 
states needed to compute the shifted energy function m 
here four major shells more than in the basis have been 
taken (e.g., in the maximum N = 2n + l of the oscilla¬ 
tor determinant is 1. Hence, for the intermediate states all 
orbits up to N=5 have been used). The results obtained 
were identical to those obtained analytically in ref. m, 
which is a good check of the convergence of the numerical 
procedure. 

In the next step then, for increasing size of the single 
particle basis up to A' = 2n + I = 18, in each nucleus 
and each basis system always three different Hartree-Fock 
calculations have been performed : 

First, a usual spherical Hartree-Fock calculation was 
done, in which the energy 

< = {D^\H\Dn) 


(25) 
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Fig. 10. Same as in Fig. 7, but for the nucleus ^®Si. 


is minimized and after convergence corrected by subtract¬ 
ing the expectation value of the kinetic energy of center 
of mass motion 

- 2 

En = E'^ - {Dn\^\Dr.) (26) 

This is the normal approach as indicated by the subscripts 
n at the total energy and the wave function. 

Second, a corrected spherical Hartree-Fock calculation 
has been done, in which the expectation value of the in¬ 
ternal hamiltonian 

is minimized. The subscript c refers to this corrected ap¬ 
proach. 

Third, a spherical Hartree-Fock calculation with pro¬ 
jection into the center of mass rest frame before the varia¬ 
tion as described in section 2 has been performed. In this 
case the energy-functional 0 

^ ^ {D^AHCmOpr) . 28 ) 

{Dpr\CmDpr) 

has been minimized. Here, for the intermediate states again 
always four major shells more than for the basis have been 
taken into account. Note, that \Dn), \Dc) and \Dpr) result 
from different variational calculations and are hence dif¬ 
ferent. 

Finally, for the largest basis system (TV = 18), we have 
studied in each nucleus a widely used approximate de¬ 
scription to deal with the center of mass motion : instead 


of minimizing E^ out of equation m one minimizes 

(* - in) + 

-I-A • {Dx\ {^2MA 

with a large Lagrangian multiplier A, i.e., one penalizes 
center of mass excitations. This prescription is exact for so 
called complete nhuj configuration spaces [2], however, is 
often applied also in truncated shell-model spaces [18,19]. 
The internal energy (i.e., ^ without the penalizing term) 
has been obtained for the three different A values 10, 100 
and 10000. 

In addition, again always for the largest basis, the 
hole-energies obtained for the corrected approach 

E^ = E,- {D,\bi{D,) (^H - 2M{A-1) ) 

(30) 

have been compared with the projected results 

j^h jp {Dpr\bj^{Dpr)HC{0)bh{Dpr)\Dpr) 

^ {Dpr\bl{Dpr)C{0)bh{Dpr)\Dpr) 

as well as the corresponding results for the simple oscilla¬ 
tor occupations. For the nuclei, in which two s-states are 
occupied (^^S and ^°Ca), obviously an additional diago- 
nalization has been performed. 

Note, that the definition (HU differs from the usual ex¬ 
pression since via the kinetic energy of the center of mass 
motion the internal hamiltonian becomes A-dependent. 
The resulting difference with respect to the usual expres¬ 
sion is for non-spurious oscillator hole-states 3?ia;/4(A—1) 
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Fig. 11. Same as in Fig. 7, but for the nucleus ®^S. 


and in the general case always larger than this lower limit. 
Inserting the results for fiuj into this formula, one obtains 
considerable effects even for the larger A-values considered 
here. 

The results for the total binding energies and the hole- 
energies of the considered nuclei are summerized in figures 
1 to 6. The left side of each figure presents the total bind¬ 
ing energy as function of the size of the basis. Three differ¬ 
ent curves are plotted : open circles refer to the Hartree- 
Fock results where the kinetic energy of the center 
of mass motion is subtracted after the variation, full tri¬ 
angles give the results of the corrected approach (1^ . in 
which this subtraction is done before the variation and 
full circles display the results of the spherical Hartree- 
Fock calculations with projection into center of mass rest 
frame before the variation 09. 

For pure oscillator occupations (i.e., the smallest basis) 
these three curves obviously coincide, for larger basis sys¬ 
tems, however, they differ considerably, i.e., display rather 
different major-shell-mixing. Let us first concentrate on 
the unprojected approaches 09 and 09 . As expected, 
the corrected approach 09 yields always a lower binding 
energy than the normal one 09 , however, for all but one 
of the considered nuclei the corresponding curves run al¬ 
most parallel with increasing basis size and their difference 
is rather small. The exception is the case of '*He, where 
09 even yields a decrease in binding energy with the ba¬ 
sis size such indicating that the underlying wave functions 
have a rather different structure than those obtained via 
09 . The energy gain of the projected approach 09 with 
respect to the corrected prescription 09 is in all con¬ 
sidered nuclei (except ^He) much larger than that of the 
latter with respect to 09 . For "^^Ca in the largest ba¬ 
sis system, e.g., the projected binding energy is 1.25 MeV 
lower than the corrected result, while the latter is only 


136 keV lower than the normal one. Note, that these 1.25 
MeV amount to almost 20 percent of the total major shell 
mixing obtained in the corrected approach 12711 . Thus, ob¬ 
viously, the restoration of Galilean invariance yields a con¬ 
siderable effect on the total binding energy and should not 
be neglected even for nuclei as heavy as "'^°Ca. 

Furthermore, as can be seen from the inverted trian¬ 
gles in the figures, the prescription 09 , which penalizes 
center of mass excitations, fails completely. For the largest 
Lagrange multiplier (A = 10000) the procedure yields in 
all considered nuclei just the simple oscillator occupa¬ 
tion. This, definitely, is a non-spurious state (i.e., contains 
no center of mass excitations), however, the major-shell- 
mixing is completely supressed in this solution. This is a 
severe warning to use the prescription 09 in incomplete 
model spaces : it always prefers non-spurious (one valence 
shell) solutions and is hence uncontrollable even if the con¬ 
figuration space is less severly truncated as in the simple 
Hartree-Fock approach discussed here. 

On the right side of figures 1 to 6 we display the proton 
and neutron hole-energies in the considered nuclei. Always 
the corrected results ll^ (indicated by the label no) are 
compared with the projected energies iImIi (indicated by 
the label pr) for both the oscillator occupation (HO) as 
well as the Hartree-Fock approach (HF) calculated in the 
= 18 basis. 

Though the underlying wave functions (and total bind¬ 
ing energies) are considerably different, in all considered 
nuclei the harmonic oscillator approach and the Hartree- 
Fock method yield remarkably similar results. For the 
non-spurious hole-states out of the last occupied major 
shell in the harmonic oscillator approach the corrected and 
projected results have to be identical as demonstrated an¬ 
alytically in ref. m and this feature holds to a large extent 
for the Hartree-Fock results, too. For the hole-states out 
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of the second and third but last occupied shell corrected 
and projected results display in both approaches rather 
similar pronounced differences. So, e.g., in the pro¬ 
jected p-holes are more than 6 MeV lower in energy than 
the corresponding corrected results and even in ^'’Ca the 
differences are still about 2.5 MeV for both the p- and the 
lowest s-holes. It was demonstrated in ref. m that these 
differences are consistent with the differences in the spec¬ 
troscopic factors out of ref. cni. This will be discussed in 
more detail in the second of the present series of papers. 

An interesting observation is made for the two nuclei 
^®Si and Here, the pl/2-holes are almost unaffected 
by the projection, while the p3/2-holes show the same dif¬ 
ferences as, e.g., observed in ™Ca. Since the coupling of 
pl/2 and d5/2 to angular momentum one is not possible 
and the d3/2-orbit is unoccupied in these two nuclei, this 
observation points to the dominance of angular momen¬ 
tum one couplings for the hole energies. 

The figures 7 to 12 demonstrate the effects on the 
charge form factors and corresponding charge densities. 
Usually, the operator for the charge density in momen¬ 
tum representation is written as CTTTI 

Pn = exp{ig • r J (32) 

r 2=1 


factor has the form 

= {D\MD) (37) 

and the corresponding charge density is just the Fourier- 
transform of this expression. 

Obviously, to obtain a translational invariant form for 
the charge density operator, (ESJ has to be complemented 
with the so-called Gartenhaus-Schwartz operator as has 
been demonstrated in ref. El 

Pinv = Pn exp{-*q • i?} (38) 

The Galilei-invariant form of the charge form factor is 
thus 

^ {D\p,nvCmD) 

{D\CmD) 

The matrix elements needed to compute this expression 
are given in appendix B. The corresponding charge den¬ 
sity is then again obtained by Fourier-transforming this 
expression. 

It has been demonstrated already some time ago [H] 
that using the so-called Gaussian-overlap-approximation 
for both the shift as well as for the Gartenhaus-Schwartz 
operator PI reduces to the “dynamically corrected” charge 
form factor 



where r is the isospin projection (proton or neutron) and 
the nucleon charge form factors fr are given by 


/.(Q2) = G^;(Q2) 


8M2 1 + 

(33) 


with the Sachs-form factors parametrized in the well- 
known dipole form (see, e.g., M) 




1 -f 




G^(Q2) = 


(843MeV)2_ 

f^p — H“2.793 
=-1.913 

1 


4M2 1-1-5 5 

1 -I- 


G^(Q2) (34) 


Here M is the nucleon mass and the negative square 
of the 4-momentum transfer 


= (ncq)^ - {AEf 


(35) 


with AE being the energy transfer and q the 3-momentum 
transfer to the system. For elastic electron scattering the 
energy transfer is given by the recoil energy {ficq )^/ {2AMc^) 
so that here 




{ficqf 


{hcqf \ 
4A2M2c4 J 


(36) 


If, as in our case, the ground state is described by a single 
determinant |D), then the “normal” elastic charge form 




Fch{Q^) exp 


3 \ 

8 {D\P^\D)] 


(40) 


In case that \D) is a non-spurious oscillator state the ex¬ 
ponential factor in (40) gets the form exp{(q6/2)^/A}, 
which is the famous Tassie-Barker correction [5]. 

On the left side of the figures 7 to 12 we compare 
for the considered nuclei the normal form factor (EH for 
the oscillator occupation (inverted open triangles) with 
the corresponding projected one II30II (full triangles), the 
dynamically corrected one p| . resulting from the solu¬ 
tion \Dc) of the minimization of the Hartree-Fock energy 
functional m, and, finally, the Galilei-invariant one PI 
computed from the solution \Dpr) of the minimization of 
the projected energy functional II28II (full circles). The cor¬ 
responding charge densities are given on the right side of 
the figures. All Hartree-Fock results have been obtained 
using the largest basis with 19 major oscillator shells. 

Let us first concentrate on the oscillator occupation. 
In "^He a large difference between the normal and the pro¬ 
jected oscillator form factors at high momentum trans¬ 
fer and consequently for the charge density at small radii 
is observed. Since we have here a non-spurious oscillator 
state this difference is entirely due to the Tassie-Barker 
correction. This correction decreases with increasing mass 
number and can in "^^Ca almost be neglected. On the other 
hand the difference of the Hartree-Fock results with re¬ 
spect to the oscillator ones increase with increasing mass 
number due to the increasing major shell mixing. E.g., in 
^°Ca Hartree-Fock and oscillator results look rather dif¬ 
ferent. 

Though computed with rather different wave functions 
the projected and dynamically corrected form factors and 












12 


R.R. Rodriguez-Guzman and K.W. Schmid: Spherical Hartree-Fock calculations with linear momentum ... 


o 

n 

E 

o 

o 

E> 

n 

£ 

U 

o 

0) 

IB 

3 

O' 

lA 



0.16 

0.14 

0.12 

0.10 

0.08 

0.06 

0.04 

0.02 

0.00 


0 1 2 3 4 5 1 2 3 4 

3 - momentum transfer * Fm r / Fm 

Fig. 12. Same as in Fig. 7, but for the nucleus ^'^Ca. 
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Fig. 13. The inelastic (mathematical) Coulomb sum rule as function of the 3-momentum transfer. Compared are the results 
of calculations using an oscillator occupation and an uncorrected elastic charge form factor (for point protons) with those using 
the same occupation but including the Tassie-Barker factor for the elastic charge form factor, those of a normal Hartree-Fock 
calculation including the dynamical correction for the elastic charge form factor and, finally, those of the projected CMPSHF 
approach. The upper part of the figure gives the results for ^He, the lower part those for the nucleus ^^0. 
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charge densities display only rather small differences in all 
the considered nuclei except "^He. This is somewhat sur¬ 
prising since in ref. laiHi larger effects of the projection 
have been seen even for nuclei up to A=40. However, be¬ 
sides being limited to projection after the variation, these 
calculations had been done with different effective inter¬ 
actions than used in the present work. Thus the present 
results do not indicate that for the elastic form factors and 
charge densities the dynamical correction is good enough 
and no projection is needed. Instead a more careful study 
using various effective interactions is required. 

In addition to the form factors we have studied the 
mathematical Coulomb sum rules, too. Details of its def¬ 
inition can be found in ref. m As usual we assumed 
point nucleons, i.e., we set the nucleon form factors out of 
eq. to 1 for the proton and 0 for the neutron. Fur¬ 
thermore, again as usual, we subtracted the square of the 
elastic form factor and divided the result by the charge 
number in order to obtain the so-called inelastic Coulomb 
sum rules. Without any center of mass correction these 
have the form 

^ | |/3„/3 J |D„) - (F",(g))^} (41) 

where the elastic form factor out of eq. has to be 
taken in the point proton limit. If we include the dynam¬ 
ical correction for the elastic form factor we obtain the 
result 

| | ^ ( 42 ) 

with the point proton limit for expression (HOI), and, fi¬ 
nally, the linear momentum projected expression has the 
form 

^inel,prr \ 4 J {I^pr\PnPn^ C {0)\Dpr) (■ cipr ^ >,2 | 

(43) 

with the point proton limit of the form factor (39). Explicit 
forms for the matrix elements entering the expressions (41- 
43) are for spherically symmetric determinants \D) given 
in appendix C. Note, that it is irrelevant whether in the 
first term of these expression the normal (32) or the in¬ 
variant form (38) of the charge density operator is used, 
since in these matrix elements the Gartenhaus-Schwartz 
operator does drop out. As for the charge form factors and 
densities the normal (41) and projected (43) results for the 
oscillator occupation have been compared with the dy¬ 
namically corrected (42) and the projected (43) Hartree- 
Fock results. For the computation of the latter obviously 
again the solutions \Dc) and \Dpr) of the corresponding 
variational calculations have been taken. 

The results for '*He and are displayed in figure 
13. Plotted are the inelastic sum rules as defined above 
as functions of the 3-momentum transfer q. As expected 
from the similarity of the dynamically corrected and the 
projected elastic Hartree-Fock form factors almost no dif¬ 
ferences between these two approaches are obtained for the 


inelastic sum rules either. That the projected oscillator re¬ 
sults almost coincide with the Hartree-Fock results, too, 
is a clear indication that the inelastic sum rule is rather 
insensitive to the major shell mixing. However, all these 
results approach the limit of 1 considerably slower than 
the normal oscillator approach. This difference is entirely 
due to the square of the elastic form factor in the above 
expressions and demonstrates that a correct treatment of 
the latter (either exact by projection or approximate by 
the dynamical correction) is definitely required. 

Since for the inelastic sum rules of the other considered 
nuclei the same behaviour as demonstrated in figure 13 is 
obtained, we shall not discuss them here. 


4 Conclusions. 

We have presented the total binding energies, hole en¬ 
ergies, form factors and charge densities as they result 
from spherical Hartree-Fock calculations with projection 
into the center of mass rest frame before the variation for 
the six nuclei “^He, ^®Si, and “^^Ca and have 

compared them to the standard Hartree-Fock results ob¬ 
tained by subtracting the kinetic energy of the center of 
mass motion either after or before the variation. Further¬ 
more, for the two nuclei '*He and we have discussed 
the inelastic Coulomb sum rule resulting from these dif¬ 
ferent approaches. For comparison, in addition the results 
for pure oscillator occupations have been discussed. 

For the total binding energies considerable effects of 
the correct treatment of Galilei-invariance are seen. In all 
the considered nuclei the energy gains of the momentum 
projected solutions with respect to the conventionally cor¬ 
rected approach using just the internal hamiltonian (which 
contains already the usual 1/A effect) in the variation are 
a considerable portion of the gains due to major shell mix¬ 
ing and hence as important as the latter. It was further¬ 
more demonstrated that the often used approximate pre¬ 
scription to penalize center of mass excitations by an ad¬ 
ditional term in the variation does not work at all at least 
in the severly truncated configuration spaces used here. 
There are strong indications that this prescription does 
only work in complete nfiw-spaces and is uncontrollable 
even if used in less severly truncated shell-model spaces. 

For the hole energies essentially the same features as 
in ref. m are observed. While the energies of the holes 
out of the last occupied shell are almost unaffected, the 
projected energies of the holes out of the second and third 
but last shell are considerably different from their conven¬ 
tionally corrected counterparts (which again include the 
trivial 1/A effects). 

For the elastic charge form factors and densities (ex¬ 
cept for the lightest considered system) there are little dif¬ 
ferences obtained between the projected and the dynam¬ 
ically corrected approach though these two approaches 
use rather different wave functions resulting from differ¬ 
ent variational calculations and the same holds for the 
inelastic Coulomb sum rules. However, these results may 
be changed, if a different (more realistic) effective interac- 
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tion is used, and hence have to be interpreted with some 
care. 

In conclusion, it has been demonstrated that a correct 
treatment of Galilei-invariance in the nuclear many-body 
problem is possible via projection methods and that its ef¬ 
fects are not only important for simple harmonic oscillator 
configurations as shown in refs. |T711llllll2l[T^ but also for 
more realistic wave functions. We shall show in the second 
of these two papers that this holds also for the spectral 
functions and spectroscopic factors. Thus we think that 
the (up to now mostly neglected) restoration of Galilei- 
invariance is unavoidable in future nuclear structure cal¬ 
culations, and, on the long run, should also be done in 
more sophisticated approaches like the shell-model [18], 
the quantum Monte-Carlo diagonalization method m or 
the VAMPIR approach |22I- 


We are grateful that the present study has been supported by 
the Deutsche Forschungsgemeinschaft via the contracts FA26/1 
and FA26/2. 


5 Appendix A: Oscillator matrix elements of 
the shift operator. 

The single particle matrix elements Sik{a) (10) of the shift 
operator 0 within oscillator single particle states play an 
essential role in the projected formalism presented in sec¬ 
tion 2. They have been given already in ref. j^]. Using the 
usual quantum numbers r, n, I, j and m for the isospin 
projection, the node number (starting from zero), the or¬ 
bital angular momentum I, which is coupled with the spin 
to the total angular momentum j and its 3-projection m, 
we obtain 


{Tinilijimi\S{a)\T 2 n 2 l 2 j 2 m 2 ) = (Jnra exp <j 


E 


i [1 + (-)'i+'=+^] 


x(_)R+/2-b]/2^(2ji + 1)(2 j2 + 

x(iij2T|l/2-l/20)(-)^'=-™^ 

X {jij2L\mi - TO2A) (a) 

(44) 


where a = \a\/b with b being the oscillator length and 


^niq„.i2(^) ^ expj + ia^l (-)"i+"2 J 

0 

XK^ jL^Ka) (45) 

where !?„;(«:) are the (dimensionless) polynomial parts of 
the usual spherical radial oscillator functions depending 
on the dimensionless variable k. An analytical form of the 
expression has been given in ref. and will not be 
repeated here. 


In case that the shift vector can be put in z-direction 
as in the spherically symmetric systems considered here, 
then _ 


47r 


2L + 1 


YUz) = 5 AO 


(46) 


Then, in eq. M . obviously, mi and m 2 have to be equal 
and the evaluation of the formulas in section 2 is simplified 
considerably. 


6 Appendix B:The projected charge form 
factor. 

In this appendix we shall give the formulas needed to 
evaluate the projected charge form factor out of eq. (39). 
Again we assume that the determinant \D) is spherical 
symmetric. This allows to fix the direction of the momen¬ 
tum transfer to the z-axis. Furthermore it can be shown 
easily, that the dependence on the angle Lpa of the shift 
vector is is trivial and can be integrated out analytically. 
Left to be done is then a two-fold integration over the 
length of the shift vector and over the angle Ha between 
the shift vector and the z-axis (direction of the momentum 
transfer). After some tedious but straightforward calcula¬ 
tion we obtain 




Jo 

2Re 


A j 

daa^ exp ^——a^ )■ / di^aSinr^a 


[ detz^(6g, a, i9)j 

T—p^n T 

E yh^h2ib<l^ a,-d) 


hi ,/i2 >0 




(47) 


where b is again the oscillator length, a = \a\/b, and h 
denotes the time reversed partner of the hole state h. The 
second sum in eq. m is restricted to positive values of 
the 3-projections of the two hole states. Furthermore q 
denotes the 3-momentum transfer while is the (neg¬ 
ative) square of the 4-momentum transfer as in section 
2. The nucleon form factors /t(( 5^) are given by eq. ijSSJ. 
Furthermore 


zUbq, a,d)=J2 V{2n + 1){2 j2 + 1 )(-)^^—^ 

L 

ijij2L\mi - TO2A)(-)'^i"i/2(jij2T|l/2 - 1/20) 


1 _ V2L+ 1 

{hhLlOOO) 2 ^ ^ ’ J {hl2L\l-lO) 


X 
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xi Tl - + ~ ^^12(^12 + 1 ) 1 

2 + 1 ) 12(12 + 1) J 

^ {LL^LilA - A0)d^lida) 

L 1 L 2 I 




(a)l 


'(2Li + 1)(2L2 + 1) 


2 I 2 + 1 


iliLil\000)^{liLl2\000){liLil\000){LLiL2\000) 


min(Zi ,L,Li) 

+2 ^ (-)^(ZiL/ 2 |A-AO) 

A^l 

ihLillX - A0)(LLiL2|A - A0)| 

(48) 


where the ry’s are given by expression (gU and 


ki 2 = (h - ji)(2ji + 1) + {h - J2)(2j2 + 1) (49) 


The matrix elements of y [2 have exactly the same form as 
except that the imaginary unit i has to be replaced 
by —i and the argument in the first rj has to be multi¬ 
plied with a factor [A — 1). Note, that the expression lli^ 
includes both natural and unnatural parity terms in the 
sum over L. The latter had been neglected in ref. [H]. 


and the 77 ’s are given by the expression 03). 

In order to evaluate the corresponding Galilei-invariant 
expression for spherically symmetric determinants \D) the 
shift vector can again be put in z-direction. We obtain 


= 


{D\PnpicmD) 

{D\CmD) 


= Z + eyipl -\{bq) 


00 

{D\C{Q)\D) [ daa^ {D\S{e, ■ a)\D) 


/ 


E 




0 

1 /^(p) 


2L-b 1 


^Acidb) PcAia) 


.AC 


I-L+1 


ABCD 


JA jc L 

JB jo I 


^ADidb)MgQ{qb) ) Pca(«)Pds(“) 


(52) 


where A, B,... denote the quantum numbers of the oscil¬ 
lator single particle basis states haIaJAi nslB^B,---, ol = 
a/b with b being the oscillator length parameter, and 


7 Appendix C:The mathematical Coulomb 
sum rule. 

In this appendix we give the explicit formulas for the ma¬ 
trix elements entering expressions for the inelastic 

Coulomb sum rules. In the normal approach one obtains 
for spherically symmetric Hartree-Fock transformations 


Sriq) = {D\pnpi\D) = Z + exp<{ -A(6g) 


([E “(2* + i)xi 

N ^h^hjh 


O^hlh', Oihl'h 


m 


- E “ E “E 


ahlhjh 


^{Ih, Ih', L){jhjh'L\l/2-l/20)^ 


{2jh + l){2jn^ + l)(^XL 


ahlhj 


(50) 


where A{lh, Ih', L) = 1 ii \ Ih - lh' \ < L < Ih + lh' and = 0 
else. 


M^ciqb) = ^ [1 + (-)'"+'^+^] A{Ia, Ic, L) 
^{2ja + l){2jc + l){-y^-^/^{jAjcL\l/2 - 1/20) 

^_'jnA+nc Ancle 

(53) 

are the reduced oscillator single particle matrix elements 
of the normal charge density operator in momentum rep¬ 
resentation. The supersript (p) at the sum symbols means 
that only proton-orbits are considered. Finally, 

PcAia) = '^{-y^~"'ijAjcL\m - mO) p’yy} (a) (54) 


where 


Xc2(a) = 

HH' 

Si^,IaSj,,ja (55) 

In eq. (E3 use has been made of the fact, that the single 
particle matrix elements do not mix different isospin pro¬ 
jections, and, for the shift vector in z-direction, do not mix 
states with different total angular momentum projections 
m either. This is indicated by the superscipts (pm). 




n n' 
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